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In this paper we demonstrate that Lindblad equations characterized by a random rate variable 
arise after tracing out a complex structured reservoir. Our results follows from a generalization of 
the Born-Markov approximation, which relies in the possibility of splitting the complex environment 
in a direct sum of sub-reservoirs, each one being able to induce by itself a Markovian system 
evolution. Strong non-Markovian effects, which microscopically originate from the entanglement 
with the different sub-reservoirs, characterize the average system decay dynamics. As an example, 
we study the anomalous irreversible behavior of a quantum tunneling system described in an effective 
two level approximation. Stretched exponential and power law decay behaviors arise from the 
interplay between the dissipative and unitary hopping dynamics. 
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I. INTRODUCTION 

The dynamic of a small quantum system interacting 
weakly with uncontrollable degrees of freedom is well un- 
derstood when a Markovian approximation applies. In 
this situation, after tracing-out the environment, the sys- 
tem density matrix evolution can be well approximated 
by a Lindblad equation pj, • 

Besides that the applicability of the Markovian approx- 
imation range over many areas of physics [E El IE IE 
there exist several real systems whose dynamics present 
strong departures from it. The main signature of this de- 
parture is the presence of strong non-exponential decay 
behaviors, such as power law and stretched exponential. 
Some examples are nanocrystal quantum dots under laser 
radiation 

SHE!, superconducting qubits [Tl 111 IT^ . 
spin environments |14|. dephasing in atomic and molec- 
ular physics 1151 . electron transfer and exciton dynamics 
in proteins |lq. and molecular systems maintained in a 
glassy environment jl7j . to name but a few. These and 
another specific experimental situations rise up the ne- 
cessity of finding formalism and effective evolutions able 
to describe the corresponding non-Markovian dynamics. 

When the environment is modeled as an infinite set 
of normal modes, departure from a Markov approxima- 
tion can be related to the corresponding spectral density 
function. This situation was extensively studied for the 
spin-boson and boson-boson models, where exacts solu- 
tions are available 0, IH l2fl l2ll |22| . Exact solutions 
can also be formulated for more general systems. Never- 
theless, due to the huge analytical and numerical efforts 
needed for getting the non-Markovian system dynamics, 
alternative numerical methods based in a special decom- 
position of the spectral density function were formulated 

Anomalous system dynamics also arise from a random 
matrix modeling of the system-environment interaction 
HE HE HE • This approach naturally applies when 
describing environments characterized by a complex dy- 
namics as for example chaotic ones. 

Outside from a microscopic point of view, there exist 



an increasing interest in describing non-Markovian effects 
in open quantum systems by introducing memory contri- 
butions in Lindblad evolutions |2£| y(J, HH • This pro- 
cedure provide easy manageable dynamics. While most 
of these models are phenomenological, in this paper we 
will relate the presence of strong memory effects in the 
standard Lindblad theory [3(j with the microscopic in- 
teraction of a system with a complex structured environ- 
ment. 

We will base our considerations in a splitting of the full 
Hilbert space of the bath in a direct sum of sub-reservoirs, 
constructed in such a way that each one guarantees the 
conditions for the applicability of a Markov approxima- 
tion. Our motivation for formulating this splitting comes 
from systems embedded in glassy environments, where 
the underlying disorder produce localized bath states, in- 
ducing a natural shell structure of modes, each set having 
a different coupling strength with the system |17|. Thus, 
we associate a different Markovian sub-bath to each set 
of states. 

As we will demonstrate, the splitting assumption al- 
lows us to generalize the usual Born-Markov approxima- 
tion, which in a natural way leads to the formulation of 
Lindblad equations characterized by a random dissipative 
rate. As is well known from a classical context, master 
equations with random rates [IE 03 EE EH are use- 
ful for describing strong non-Markovian effects [IE EU . 
Here we will demonstrate that the same scheme can also 
be applied in a quantum context under the previous con- 
ditions, i.e., a complex environment under the splitting 
condition. 

As an example, we will study the anomalous dissipative 
dynamics of a quantum tunneling system described in a 
two level approximation 0,0]. Strong non-exponential 
behaviors, such as stretched exponential and power law, 
arise from the interplay between the unitary hopping dy- 
namics and the memory effects induced by the environ- 
ment. The conditions under which our modeling can be 
mapped to the spin-boson model and stochastic dynam- 
ics are established. In this context, the differences be- 
tween our framework and other approaches p3l l24j in- 
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troduced to deal with non-Markovian environments are 
established. 



II. REDUCED DYNAMICS FROM COMPLEX 
ENVIRONMENTS 



In general, the evolution of a system interacting with 
a complex environment can not be described in a Marko- 
vian approximation. While a general treatment is clearly 
not possible, when the splitting condition applies, for 
weak system-bath interactions, the reduced dynamics 
can be described through a generalization of the Born- 
Markov approximation. 



A. Generalized Born-Markov approximation 

We start by assuming a full microscopic Hamiltonian 
description 



Ht — Hq 



(1) 



where Hs correspond to the Hamiltonian of a system S, 
and Hb correspond to the Hamiltonian of the bath B. 
The term Hj = qs ® Qb describes their mutual interac- 
tion, with the operators qs and Qb acting on the system 
and bath Hilbert spaces respectively. 

In an interaction representation with respect to H$ + 
Hb, the total density matrix pr(t) evolves as 



dprjt) -i, . . , 
-^ r = T [H I (t),p T (t)], 



(2) 



where Hi (t) is the interaction Hamiltonian in the Hciscn- 
berg representation. Integrating formally this equation, 
and substituting the solution in it, the evolution of the 
reduced system density matrix ps(t) — Trs{pT(t)} can 
be written as 



dps® 
dt 



dt'Ti B {[Hi(t),[Hi(t'),p T (t')}]}, 

(3) 

where, as usual, a first order contribution was discarded 
after assuming Tr B [Hi(t)pT{0)] = 0. From this evolu- 
tion, the well known Born-Markov approximation can be 
deduced 0. Q. The Born approximation consists 
into assume, at all times, an uncorrelated structure for 
the total density matrix 



Pr(t) « p s (t) <g> p B , 



(4) 



where ps define a stationary state of the bath. This 
assumption is consistent up to second order in the inter- 
action Hamiltonian. When the decay of the bath corre- 
lation defines the small time-scale of the problem, after 
introducing Eq. in Eq. ©, the Markovian approx- 
imation leads to a closed local in time density matrix 
evolution. 



We remark that the Born-Markov approximation does 
not rely in any specific model of bath dynamics |ES>IEE1> 
such as an infinite set of harmonic oscillators. In fact, 
its master equation is independent of model assumptions 
used in its derivation 0. 

Now we consider a complex environment for which the 
previous approximations are not valid. As is usual when 
dealing with complex environments |25l l26l 1271 128( , in- 
stead of defining the bath Hamiltonian Hb as an infi- 
nite set of normal modes, here we specify it through its 
eigenstates basis {|e)}, which in a weak interaction limit, 
results unmodified by the interaction with the system. 
As a central hypothesis, we will assume that, while the 
full action of the environment can not be described in 
a Markov approximation, it is possible to split the full 
Hilbert space of the bath as a direct sum of subspaces, 
in such a way that each one defines a sub-reservoir able 
to induce by itself a Markovian system evolution [3f^ |. 
These hypotheses are the main assumptions that allow 
us to formulate our results. 

In conformity with the splitting condition, we write 
the interaction Hamiltonian as a direct sum of sub- 
Hamiltonians 



Ht = Ht 



Ht 



Hi 



(5) 



where Hi R = qs <8> Qb r - Here, each operator Qb r de- 
fines the interaction between the system and each sub- 
reservoir R. 

In order to describe the joint action of all sub-reservoirs 
over the system, instead of the uncorrelated form Eq. 10}, 
we introduce the generalized Born- approximation 



p T (t)ttJ2p R (t)®E R , 



where Tr s[PR(t)] = 1, and we have defined 
Zr = ( e «l Pb \*r) 

{tR} 



(6) 



{\£r)} is the base of eigenvectors that span the subspace 
corresponding to each sub-reservoir. Therefore, each con- 
tribution in Eq. JBJ consists in the external product be- 
tween a system state p R (t) and the projection of the sta- 
tionary bath state ps over each subspace R. In physical 
terms, each state p R {t) takes in account the dissipative 
effects induced by each sub-reservoir. 

After introducing Eq. © in Eq. we get the ap- 
proximated evolution 



dps® 
dt 



dt' (8) 



Tr BR {[Hi R (t), {Hi R (t'),p R (t')®p BR }}}, 



where Trs R {»} means a trace operation with the states 
{|£-r}} corresponding to each subspace. Furthermore, 
we have introduced the sub-bath density states Pb r — 
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Sr/Pr, where 



Pr = Ti Br {^r} = ( e n\ Pb \zr) 



(9) 



The normalization condition Tr b[pb] = 1 implies the 
relation J2rPr = 1- Thus, the set {Pr} can be seen 
as a set of probabilities defined by the weight of each 
subspace in the total stationary bath state. 
From Eq. © we can write 



Ps (t) = Tr B [p T (t)}ttJ2PRpR(t)- 



(10) 



Then, the evolution Eq. JHJ is in fact a linear combination 
of the evolutions corresponding to the set {pR(t)} : each 
one participating with weight Pr. After introducing the 
Markovian approximation [EE IE El to the evolution of 
each state PR(t), in a Schrodinger representation, we get 



dpR(t) -i 

~aT = T [Hs > PR{t)] 



dt' 



h / JO 

Tt Br {[H Ir: [H lB {-t'),p R {t)®pB H ] 



(11) 



This evolution correspond to the usual Born-Markov ap- 
proximation when considering a bath consisting only of 
the subset of states {|cr)} and characterized by the sta- 
tionary state pb r ■ The system density matrix is defined 
by the linear combination Eq. I|10|) . 



B. Random Lindblad equations 

The evolution Eq. I|ll|l. disregarding transients of the 
order of the sub-bath Hamiltonian correlation time, can 
be always well approximated by a Lindblad equation 



dpRjt) 
dt 



C H [p R (t)}+ lR C[p R (t)}, 



(12) 



where = (—i/H)[Hs,»] is the system Liouville su- 

peroperator and the Lindblad superoperator is defined 
by 



(13) 



As the underlying microscopic interaction between the 
system and the environment is the same in each sub- 
space, the set of operators {V a } does not depend on index 
R. Nevertheless, each subspace has associated a different 
characteristic dissipative rate j R . As this rate arises from 
the interaction of the system with the manifold of states 
{|cr)}, consistently with the Fermi golden rule Q, it is 
proportional to the characteristic interaction strength of 
each subspace, denoted as \Qb r \, multiplied by the corre- 
sponding density of states g R (e) = J2{e R } ^( e ~ e n) eval- 
uated in a characteristic frequency ujs of the system, i.e., 

7R ~ \QB R \ 2 9R{f^Js)- 



With these definitions in hand, we conclude that un- 
der the generalized Born-Markov approximation, we can 
represent the dynamics induced by the complex environ- 
ment by a Lindblad master equation characterized by a 
random rate variable, defined by the set {7r,Pr}. Cor- 
respondingly, the system state follows from the average 



p S {t) = Y,PRPR{t) = <Pfl(*)>. 



(14) 



Random rate equations were extensively used to model 
classical anomalous diffusion processes in disordered me- 
dia [IE HE HEl- Here, we have derived a similar 
structure for a different physical situation, i.e., quantum 
systems embedded in a complex structured environment. 

As in a classical context, while the density matrixes 
PR{t) follows a Markovian evolution, the average state 
ps(t) evolve with a non-Markovian evolution. This evo- 
lution can be easily obtained in a Laplace domain, where 
the average Eq. (|14fl takes the form 



Ps(u) 



1 



u - (C H + IrC) 



ps(0) = (G R (u))ps(0), 



(15) 

with u being the Laplace variable, and we have used 
the solutions p R (t) = exp[(£# + 7r£)£]ps(0). Consis- 
tently with the uncorrelated initial condition pr(0) = 
Ps(0) (g) ps(0), there not exist any statistical correla- 
tion between ps(0) and the random rate set. Thus, 
the average evolution can be obtained without appeal- 
ing to a projector technique [3 0, HE In fact, 
after introducing in Eq. I|15|) the identity in the form 
(G R (u)[u — (C R + 7r£)]) _1 , it is immediate to get the 
deterministic, closed, non-Markovian evolution equation 



dpsjt) 
dt 



= C H [p s (t)}+ f drh(t-T)[ Ps (T)}, 
Jo 



(16) 



where the superoperator L is defined in the Laplace do- 
main by the equation 



(G R (u) lR C)[*} = (G R ( U ))L(u)[. 



(17) 



Depending on the set {7r,Pr}, Eq. l(TB|l may lead to 
the presence of strong non-Markovian decay behaviors in 
the system dynamics. This characteristic originates from 
the entanglement of the system which each sub-reservoir, 
situation explicitly introduced by Eq. (JSJl. 

An example of complex structured environment where 
the generalized Born-Markov applies straightforwardly is 
a bath Hamiltonian whose eigenvectors can be labelled 
with two indexes (E, R). The index E is continuous, and 
for each R the corresponding sub-manifold of states is 
able to induce a different system Markovian-decay. The 
difference between the Markovian dynamics may origi- 
nates in the coupling strength of each sub-manifold with 
the system. On the other hand, it may originates due 
to strong variations of the bath density of states with 
index R. The system dynamics follows as a superpo- 
sition of Markovian dynamics whose weights are taken 
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in account through the generalized Born approximation 
Eq. |JBJ. If the decay induced by each sub-manifold is the 
same, the generalized Born approximation reduces to the 
usual one, and then a Markovian evolution is obtained. 
Further examples can be established in the context of 
random band-matrix models ,7|, where the Markovian 
sub-baths, for example, may be associated to subspaces 
with a different characteristic bandwidth. 



C. Effective approximation 

Classical master equations with random rates are char- 
acterized by equations similar to those obtained previ- 
ously. Nevertheless, as in general the underlying numbers 
of states is infinite, some kind of approximation is neces- 
sary in order to obtain the operator L, as for example an 
effective medium approximation |35ll36|. Here, we intro- 
duce a similar approximation in order get a general char- 
acterization of the dynamics. Thus, in Eq. H7fl we dis- 
card the dependence introduced by the Lindblad super- 
operator C in the propagator Gr(u), i.e., C — > —I. Then, 
we get the approximated solution L(-u) ~ K(u — £#)£, 
from where it follows the evolution 



d Ps (t) 
dt 



C H [ps(t)}+ / dTK(t~r)e^ c "C[p s (T)}. 
Jo 

In this approximation all information about the random 
rate is introduced through the kernel function 



K(u) 



lR 



1R 



1 



u + j R 



(19) 



As in a classical context, this kernel can be associated 
with a waiting time distribution w(t) and a survival prob- 
ability Po(t) defined by 



w(u) 



1R 



u + 7r 



Po(u) 



1 



u + 7_r 



(20) 



In classical master equations, these objects define a con- 
tinuous time random walk j^l El El El El El In 
the quantum case, a similar stochastic dynamics can 
be constructed [30j . It consists in the application at 
random times of the superoperator £ = £+1, implying 
the transformation p — > £[p], while during the inter- 
vals between these disruptive actions the system evolves 
with its unitary dynamics, U(t) — exp[i£#]. The in- 
tervals between the successive applications of £ follows 
from the waiting time distribution w(t). The func- 
tion Pq (t) defines the corresponding survival probability, 
P (t) = 1 - Jo drw{T). Thus, the average over different 
realizations of the random times can be written as 

p s (t) = P (t)e tc " Ps (0) 

+ f dTw(t-T)e^ CH £[p s (T)}. (21) 
Jo 



From here, in a Laplace domain, it is straightforward 
to recuperate the evolution Eq. l(T%|) . When C ^ £—1, 
with £ a completely positive superoperator Q , a similar 
stochastic dynamics can be formulated after introducing 
a limit procedure |40l |. 

We remark that the stochastic interpretation [Eq. I|21|l ] 
was constructed after associating to the kernel K(t) 
a waiting time distribution and a survival probability, 
Eq. if^UI) . This association does not rely in the gener- 
alized Born-Markov approximation, neither it was de- 
duced from a conditional continuous time measurement 
theory Therefore, it is not clear if one can associate 
to the stochastic dynamics a random signal of a mea- 
surement apparatus. If this is the case, contradictions 
between environmental decoherence and wave-function 
collapse may arise • 



III. QUANTUM TUNNELING SYSTEM 
DRIVEN BY A COMPLEX ENVIRONMENT 

As an example of our formalism, in this section we 
will characterize the dissipative dynamics of a quantum 
tunneling system described in a two level approximation 
|l8ll2]| and driven by a complex environment. Then, the 
system Hamiltonian can be written as 



hA 



(22) 



The first term, proportional to the z-Pauli matrix a z de- 
fine the energy of the effective levels, and the second one, 
proportional to the x-Pauli matrix <r x , introduce the re- 
versible hopping between the two effective states. 

The complex environment will be represented by the 
Lindblad superoperator 

£[»} = ^([<T z ;a z } + [a z ,»<T z ]), (23) 

and an arbitrary set {7^, Pr} of random rates and 
weights. For fixed rate, this superoperator induces a dy- 
namics equivalent to a thermal environment in a high 
temperature limit |2^ . 

The evolution of the system density matrix is defined 
by Eqs. I|ltj|) and l|17|l . Here, we write the evolution in 
terms of the components of the Bloch vector, which are 
defined by the mean value of the Pauli matrixes, Sj(t) = 
Trslps (t)<jj}, with j = x, y, and z. We get 



' / ' S ' x(/| = -w A S Y (t)- J' ' dr{T x (t-r)S x (r) (24a) 
-T(t-r)S Y (r)}, 

= tu A Sx(t) - AS z (t) - [ dr (24b) 
Jo 



dt 

dS Y (t) 
dt 

dSzjt) 
dt 



{T Y (t - t)Sy(t) + T(t - t)S x (t)}, 
AS Y (t). (24c) 
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Thus, the system evolution is completely characterized 
by three memory kernels Tx(t), IV (i), and T(t). In Ap- 
pendix A, we give the exact expressions of these kernels 
for arbitrary random rates, joint with the kernels that 
arise from the effective approximation Eq. (JTSJ. From 
Eq. (|24J) it is straightforward to write the system density 
matrix evolution [Eq. i|16|l ] as a sum of Lindblad super- 
operators, each one characterized by a different kernel. 



random rates and their corresponding weights. With 
these definitions, it is simple to demonstrate that after a 
transient of order 70, the waiting time distribution and 
its associated survival probability, Eq. il2U|) . present a 
power law decay behavior |43j, w(t) ~ l/(7oi) 1+Q , and 
Po(t) « l/(7oi) Q , where a — a/b. Clearly, this behavior 
is reflected in the system dynamics. 



A. Dispersive limit 

When the hopping frequency is zero, A = 0, the dy- 
namics reduce to a dispersive one. Thus, the coherences 
decay continuously while the population of each effective 
level remains constant. In this limit, from Appendix A, 
for arbitrary set {jr,Pr} we get the exact kernels 



r*(t) 

!>(*) 
T(t) 



K(t) cos[usAt], 
K(t) cos[ujAt], 
K(t) sm[uiAt]. 



(25a) 
(25b) 
(25c) 



where K(t) is defined in the Laplace domain by Eq. 1)190. 
We note that these kernels also arise from the effective 
approximation Eq. (|18J) . indicating that for A = 0, both 
evolutions coincide. 

From Eqs. I|24|) . the exact solution of the Bloch vector 
is given by 

Sx(t) = P {t){cos[w A i\Sx(0)-sm[uAt]S Y (0)}, (26a) 
S Y (t) = P Q (t){sm[Lo A t]S x (0) +cos[wA*]5y(0)}, (26b) 
S z (t) = S z (0), (26c) 

where Po{t) is the survival probability defined by its 
Laplace transform Eq. (|20|) . which in the time domain 
reads Po{t) = J2 R Pr exp[— 7r<]. Consistently, we note 
that the exact solutions Eqs. (|26|l correspond to an aver- 
age over Markovian solutions, each one characterized by a 
rate 7^ and participating with weight Pr. Depending on 
the distribution of the dissipation rate, arbitrary forms 
of the decay can be obtained from this average over ex- 
ponential functions. Hence the non-Markovian behavior 
can be observed in the relaxation of the density matrix 
to the stationary state. 



B. Anomalous decay behaviors 

The form of the set {7^, Pr} depends on the specific 
structure of the complex environment. Here, we will de- 
termine this set in a phcnomcnological way as a function 
of the system decay behavior. We will be interested in 
obtaining anomalous decay dynamics such as power law. 
A possible set consistent with this decay is 

7K = 7 exp[-&R], Pfl = (l-e- Q )exp[-ai?], (27) 

where R G [0, 00], 70 scale the random rates, and the 
constants b and a measure the exponential decay of the 



a q0 _ I I I I 1 1 ■ | I ■ I !■■■■[ 1 1 I I I I ll| 1 1 I I I I ll| 1 1 I I I I M 



r° 10 - 




FIG. 1: Survival probability. From top to bottom, the pa- 
rameters are /3/-y = 0.75, 10 _1 , 10 -2 , 10 -3 , 10" 4 , and the 
Markovian limit /3 = 0. In all cases we take a = 1/2. 

When < a < 1, the kernel K{t) corresponding to the 
set Eq. (|27J) can be well approximated by the expression 



K{u) 
with the definitions 

7 = (iR), 



7 



1 + 03/u) 1 -"' 



ill) 



(28) 



(29) 



The scaling of these parameters can be motivated by con- 
sidering a two dimensional set of random rates |4^ | . From 
Eq. (|19|l and l|2U|l . the waiting time distribution and its 
associated survival probability can be obtained as 



w(u) 



K(u) 
+ K(uY 



Po(u) 



K{u) 



(30) 



From here, is it is simple to proof that w(u) is a com- 
pletely monotone function [3(|, which implies that Po(t) 
decays in a monotonous way or equivalently, w(t) > 0. 

In Figure 1 we plot the survival probability Po(t) by 
assuming the kernel Eq. I|28(l for different values of (3/j. 
We note that in a short time regime, the decay is an 
exponential one, while in an asymptotic regime a power 
law behavior is present 



P (t)^exp[- 7 t], 



l-a 



1 



-yT(l - a) t c 



(31) 
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where T(x) is the gamma function. These asymptotic 
behaviors follows immediately from Eq. (|3L)[) . When 
the dispersion of the random rate 7^ is zero ((3 = 0), 
consistently the dynamics reduce to a Markovian one, 
K{u) = 7, which implies the pure exponential decay 
Po(t) = cxp[— 7*] and w(t) = 7exp[— 7*]. 

In the next subsection we will characterize the tunnel- 
ing dynamics by assuming a complex environment char- 
acterized by the random rate set Eq. I|27[) or equivalently 
by the kernel Eq. (|2"gj) . 



C. Tunneling dynamics 

Here we will analyze the tunneling dynamics for a sym- 
metric case uja = 0, which arise when the two effective 
levels have the same energy. From Appendix A, the exact 
kernels read 



T Y (u) 
T(u) 



K{u), 

K(u + A 2 /u), 
0. 



(32a) 
(32b) 
(32c) 



As before, the kernel K(u) is defined by Eq. (|19|l . The 
exact solution of the Bloch vector can be obtained in a 
Laplace domain. We get 



S x {u) = 

S Y {u) = 
S z (u) = 



1 



:Sx(0), 



1 , , • > "/• (33a) 

A(u){uS Y (0) - AS z (0)}, (33b) 
A(u){[« + T Y (u)]Sz(0) + AS Y (0)}, (33c) 



where we have defined 



A(«) 



uTy{u) + A 2 ' 



(34) 



which can also be expressed as A(u) = Pq(u + A 2 /u)/u. 

In this case it is not possible to find in the time domain 
a general exact solution for arbitrary memory kernels. A 
simple analytical solution is only available in a Markovian 

case [K(u) = 7] 



Sx(t) = 
S Y (t) = 

Sz(t) = 



e -7t/2{5 y ( ) CO sh[At] 



(35a) 



A- i [( 7 /2)5y(0) + AS z (0)]sinh[At], (35b) 
e^ t/2 {S z {0)cofih[Xt} 

+\- 1 [( 1 /2)S z (0) + AS Y (0)} sinhfAt], (35c) 



where A = \J (7/2) 2 — A 2 , and 7 defines the unique dissi- 
pative rate. Notice that in the limit of null dissipation, a 
periodic hopping between the effective levels is obtained. 

For arbitrary random rates {7.R, Pr}, the dynamics can 
be characterized in different regimes. First, in the case 



N 0.0 




At 



FIG. 2: Average Sz(t) considering as initial condition the 
upper eigenstate of a z . The envelopes are given by ±Po(t/2). 
From top to bottom, the parameters are 7/ A = 0.05, 0.15, 
and 1.0. In all cases we take 01 = 1/2, j3 = 7/2, and uja = 0. 



A » (7a), from Eqs. it is possible to get the ap- 
proximated solutions 



Sx(t) = P (t)S x (0), 

S Y {t) ~ P (t/2){cos[At]S Y (0) 

S z (t) ~ P (t/2){sm[At]S Y (0)- 



(36a) 

■ sia[At]S z {0)}, (36b) 
cos[A£]S z (0)}. (36c) 



Thus, the dynamics consist in a periodic tunneling be- 
tween the two effective states, and whose decay can be 
written in terms of the survival probability. As in the pre- 
vious case, this solution correspond to an average over the 
corresponding Markovian solutions, i.e., Eq. I|35|l written 
in the limit of small decay rate when compared to the 
tunneling frequency A. 

In Figure 2 we plot the average of the z-Pauli matrix 
which follows from Eq. with the kernel Eq. (gSJl . As 
initial condition we take the upper eigenstate of a z . We 
verified that the exact solutions are well described by the 
approximation Eq. (|36|l for parameters values satisfying 
7/A < 1. As the envelope decay is given by Po(t/2), by 
increasing the average rate 7, the dynamics decay in a 
faster way. This dependence is broken when the average 
rate is much greater than the hopping frequency. 

In the limit A <§; {jr), the dissipative dynamics dom- 
inates over the tunneling one. In Figure 3 we plot Sz(t) 
[Eq. H33|) ] for different values of the characteristic param- 
eters of the kernel Eq. (|28|l . We note that by increasing 
the average rate 7, a slower decay is obtained. Thus, the 
dynamics develops a Zeno-like effect H^, . From the 
exact solution Eq. (JJSJl, the characteristic decay of the 
Bloch vector can be approximated by the expressions 

S x (u) = P (u)Sx(0), (37a) 
S Y {u) ~ Z(u){uS Y (0) - AS z (0)}d(u), (37b) 
S z (u) ^ Z(u){S z (0)+AS Y (0)d(u)}, (37c) 
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where we have introduced Z(u) = u 1 w(A 2 /u) and the 



function d(u) = [u + K(A 2 /u)]~ 



1 P (A 2 /u). For 



the kernel defined by Eq. 1)28(1 . the characteristic decay 
Z(u) results 



Z{u) 



1 



u + Ci + C a u l 



I — a A 2ct 
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(38) 



As can be seen in Figure 3 (dotted line) , besides the oscil- 
latory behavior, after the transient jt <C 1, this function 
provides an excellent fitting of the decay dynamics. 




~i 1 1 1 1 1 r 

100 150 200 250 300 

At 



FIG. 3: Average Sz(t) considering as initial condition the 
upper eigenstate of cr z . The fitting decay curves (dotted lines) 
are given by Eq. 1381 . From top to bottom, the parameters 
are 7/A = 200, 100, 50, 25, 10, and 2.5. In all cases we take 
a = 1/2, j3 — 7/2, and ua = 0. In the inset we show the same 
graphic in a log- log scale. 

The function Z{t) is characterized by a reach vari- 
ety of behaviors. First, we note that in the Markovian 
limit, (3 = 0, we get an exponential decay with rate 
C\ = A 2 /7, which clearly diminish by increasing 7. In 
the non-Markovian case, in a short time regime, we can 
approximate 



Z(t) ~ exp{-[dt 



C a t a 



■a 



}}; 



while in an asymptotic long time limit we get 
(l-a)C a 1 



Z(t) 



r(a) c 2 t 2 -*' 



(39) 



(40) 



Thus, the dispersion of the random rate (measured by fi) 
induce, at short times, an extra stretched exponential de- 
cay, while in the asymptotic regime it scales a power law 
behavior [C 2 /C* Q = A^-^/T/ 91- "]- The characteristic 
rates of both regimes arise from a competence between 
the unitary and dissipative dynamics. We notice that by 



increasing the dispersion rate /3, the characteristic rate 
of the stretched exponential decay is increased, while the 
rate for the power law regime is decreased. The depen- 
dence in the hopping frequency A is the inverse one. 

The Zeno-like effect can be qualitatively understood 
in terms of the stochastic evolution corresponding to 
the effective approximation Eq. I|18l) . This stochastic 
process develops in the system Hilbert space and con- 
sists in the application at random times of the super- 
operator £ = C+l, which in view of Eq. i|23|) reads 
£[•] = a z ua z , while in the intermediates times the system 
evolves with its unitary evolution U (t) = exp[— iAta x /2]. 
The superoperator £ implies the disruptive transforma- 
tions Sx — * — Sx, Sy — ► — Sy, Sz — > Sz, while the 
unitary dynamics is equivalent to a rotation around the 
x-direction. In the limit of vanishing hopping frequency 
A, the continuous applications of the superoperator £ kill 
the x-y components and frozen the dynamics in the initial 
condition Sx(0). Thus, a pure Zeno effect is recuperated. 
For A/(7#) <C 1, the decay dynamics is determined from 
the competence between the transformations induced by 
£ and U(t), defining the Zeno-like regime. This interpre- 
tation is exact in a Markovian limit and always valid for 
the effective master equation Eq. (I18|) . 



D. Anomalous decay behavior from a finite set of 
random rates 



In obtaining the previous results we have assumed an 
infinite set of random rates, Eq. I|27|). whose effects can 
be approximated by the kernel Eq. (|28|1 . While this elec- 
tion guarantees the presence of an asymptotic power law 
decay, strong non-exponential behaviors can be obtained 
in an intermediate regime by considering only a finite set, 
1 < R < A r max , of random rates {7^, Pr}. On the other 
hand, for a finite set, the asymptotic system dynamics is 
always Markovian and characterized by the inverse rate 
{1/jr). This result follows from \\m u ^ K(u) — (1/jr). 

In Figure 4 we show the decay dynamics induced by 
an environment characterized by a finite set of random 
rates j R (iV max = 7) with equal weights, P R = l/N max . 
Each curve follows from a superposition of Markovian 
solutions, Eq. jSSjl with 7 — > 7^. The set of rates 
{ir} of each plot differ in a multiplicative factor, in 
such a way that the relation between the average rate 
7 = (7r) and the corresponding fluctuation rate (3 = 
\{'i R ) — {ir) 2 \/ (7r) remains constant in all curves. For 
the case 7/A = 2.5, the random rates are Jr/A = 0.59, 
1.0, 1.09, 1.21, 4.0, 4.7 and 4.88. 

In order to enlighten the intermediate non-exponential 
regime, we have plotted the shifted average S^(t) — 
a + (1 — a)Sz{t), with a <C 1. In the deep Zeno-like 
regime [7/A > 10], S^{t) can be well approximated 
by an stretched exponential behavior 5f (t) ~ a + (1 — 
a)exp[-((t) s ], with 6 w 0.7 and (/A S (0.025,0.12). 
For 7/A < 10, a power law fitting is more adequate 
S a z {t) ~ a + (1 - ct)[1 + (t}- s , with 5 6 (3.5,4.5) 
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FIG. 4: Average Sf (t) = a + (1 - a)Sz(t), considering as 
initial condition the upper eigenstate of a z and a finite set 
of random rates with equal weights. From top to bottom, 
the parameters are 7/A = 50, 25, 10, 5, 3.5, and 2.5. In 
all cases we take a — 0.01, (8/7 = 0.51, and loa = 0. For 
7/A > 10, the dotted lines correspond to the fitting S%(t) ~ 
a + (1 — a) exp[— (Ct) 5 ], while for 7/A < 10, they corresponds 
to S z (t) ~ a + (1 - cr)[l + (t]~ s (see text). 

and (/A £ (0.05,0.125). We note that a similar non- 
exponential fitting was found in Ref. 17] by considering 
the action of a finite bath, which can be associated with 
a glassy environment. On the other hand, the oscillatory 
effects in the decay of S z (t) arise from the Markovian 
solutions, Eq. I|35[) . corresponding to the rates satisfying 
7fi/A < 2. In fact, for the Markovian solution, this con- 
dition delimits the change between a monotonous and 
an oscillatory decay behavior. Consistently, we notice 
that by increasing the average rate, the amplitude of the 
oscillations are smaller. A similar effect can be seen in 
Figure 3. 

E. Mapping with other models 

Our formalism relies on the applicability of the gener- 
alized Born-Markov approximation. Here we explore the 
possibility of mapping its dynamics with other models 
that also induce anomalous decay behaviors. 

Spin-Boson model: The spin-boson model is defined 
by the total Hamiltonian 

H T = ^{uj A <J z + Aa x } + ^a z Q B + H B , (41) 

where the bath Hamiltonian Hb = YlAipj /^ m j) + 
mjUjjqj] corresponds to a set of harmonic oscillators, and 
Qb = K j1j- The bath is characterized by the spec- 
tral density function 

J(w) = (7r/2)d£^/"WW«> ~ *>i), (42) 



and assumed to be in equilibrium at temperature T. As 
is well known, the reduced system dynamics can be ob- 
tained in an exact way [H, [13, 120, ull |22| . It reads 

Sx(t) = [ dT[Yj, s) (t-T) + Y\ a \t~T)S z (T)} (43a) 
Jo 

+Y B s \t)S x (0) + Y B a) (t)S Y (0), 
SA, = (-) 

^ = j\T[Kf(t-r)-K${t-T)S z {T)] 

+K ( B a \t)S x (0)+K B s \t)S Y (0)}, (43c) 

where the corresponding kernels can be written as func- 
tions of J(w). On the other hand, it is possible to 
write the exact averaged evolution Eq. 1241) in the form 
Eqs. I|43|l . In Appendix B we present the kernels corre- 
sponding to each dynamics. From these expressions, it is 
simple to proof that to first order in A, after disregarding 
a phase factor, both set of kernels can be mapped under 
the condition 

£ P R exp[- 7fl t] = exp[-Q'(t)], (44) 

R 

where 

Q ( t ) = J-J o dw^-coth(—)[l-co S (wt)}, (45) 

define the real part of the thermal bath correlation. We 
remark that the mapping Eq. I|44|l is only valid in a high 
temperature limit, condition consistent with the Lind- 
blad structure Eq. J22J. 

In this context, from Eq. (|44Jl . it is possible to enlighten 
the difference between the present approach and that de- 
veloped in Refs. [23l[24|. In our approach, which relies in 
splitting the Hilbert space of the bath as a direct sum of 
subspaces, exp[— Q'(t)] is written as a sum of exponential 
functions, each one associated to each Markovian sub- 
reservoir. Instead, in Refs. [Ul^, Q'(t) is expressed 
as a sum of exponential functions. This representation 
relies in an artificial discomposing of the spectral density 
function J(u>) as a sum of individuals terms. Thus, the 
Hilbert space of the bath is effectively split in an exter- 
nal product of subspaces, each one associated to a non- 
Markovian sub-reservoir. As in our approach, the system 
density matrix can be written in terms of a set of auxil- 
iary sub-density-matrixes. Nevertheless, their evolution 
involves coupling among them all. 

Stochastic Hamiltonian: Decoherence in small quan- 
tum systems is also modeled by introducing stochastic el- 
ements in the system evolution. This situation arises nat- 
urally in many physical systems [Til IT2I Il3l W\ . Consis- 
tently with the spin-boson model we consider the stochas- 
tic Hamiltonian 

Hst{t) = |{k4 + £(t)K + Act,}, (46) 
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where £(t) is a classical non- white noise term. 

By assuming ((£(£)}}£ = 0, where ((•••))$ means an 
average over realizations of the noise, in the limit of van- 
ishing A it is simple to solve the stochastic dynamics 
and obtain the average of the Pauli matrixes. The fi- 
nal evolution is the same as in Eq. I|2t)|l after replac- 
ing Po(t) with the average dephasing factor D(t) — 
((exp[i J dr^r)]))^. Thus, the generalized Born-Markov 
approximation can be mapped to the stochastic Hamilto- 
nian evolution under the condition D(t) = Po(t), which 
explicitly reads 



exp[i / cLt!;(t)] 
o 



(47) 



This condition can be consistently satisfied if the dephas- 
ing factor D(t) decays in a monotonous way. 



IV. SUMMARY AND CONCLUSIONS 

We have presented a theoretical approach intended 
to describe the dynamic of small quantum systems in- 
teracting with a complex structured environment. Our 
formalism is based in an extension of the well known 
Born-Markov approximation, which relies in the possi- 
bility of splitting the environment as a direct sum of 
sub-reservoirs, each one being able to induce by itself 
a Markovian system dynamics. Then, we have demon- 
strated that the full action of the complex environment 
can be described through a random Lindblad master 
equation. The set of random rates follows from a Fermi 
golden rule. Thus, they are proportional to the char- 
acteristic coupling strength of each subspace multiplied 
by the corresponding sub-density of states evaluated in 
a characteristic frequency of the system. The associated 
probabilities are defined by the weight of each subspace 
in the stationary state of the bath. 

From a phenomenological point of view, the set of ran- 
dom rates and weights can be determined in a consistent 
way in function of the system decay. In fact, the sys- 
tem dynamics is characterized by a non-Markovian mas- 
ter equation that in function of the random rate set can 
develop strong non-exponential decays. 

As an example we worked out the dissipative dynamic 
of a quantum tunneling system in a two level approxima- 
tion. We have introduced a set of random rates that lead 
to the presence of asymptotic power law decay. In the 
limit of small hopping frequency, when compared with 
the average rate, we have showed that a Zeno-like phe- 
nomenon arises, which is characterized by a stretched ex- 
ponential and a power law decay. These behaviors follow 
from the interplay between the unitary dynamics and the 
entanglement-memory-effects induced by the reservoir. 

For the tunneling dynamics, we have also demon- 
strated that non-exponential decays arise even by con- 
sidering a small set of random rates. Furthermore, we 
have established the conditions under which the random 



Lindblad evolution can be mapped to a spin-boson model 
and a stochastic Hamiltonian evolution. 

Finally, we want to emphasize that the present results 
define a new framework for describing anomalous quan- 
tum system dynamics, which consists in taking the char- 
acteristic rate of a Lindblad equation as a random dis- 
tributed variable. We remark that this approach was not 
derived from an ensemble of identical systems whose local 
interactions with the environment can be approximated 
by different Markovian evolutions. In fact, the underly- 
ing microscopic physics can be related to a single quan- 
tum system coupled to an environment with a complex 
structured spectral density function and whose dynami- 
cal influence over the system can be approximated by a 
direct sum of Markovian sub-reservoirs. Thus, our ap- 
proach may be relevant for the description of anomalous 
decay processes in individual mesoscopic systems embed- 
ded in a condensed phase environment 0, 0, 0] . A nat- 
ural example for which the generalized Born-Markov ap- 
proximation may applies are glassy reservoirs, where the 
underlying configurational disorder produce a hierarchi- 
cal distribution of coupling strength between the single 
system and the corresponding localized eigenstates of the 
reservoir |17| . 



APPENDIX A: EXACT KERNELS 

Here we present the exact expressions for the kernels 
Tx(u), Ty(u), and T(u) that define the evolution of the 
Pauli operators average, Eqs. I|24|) . For arbitrary rates 
{"fR: Pr}, the kernels read 

T x (u) = D{[u(u + C) + A 2 }(u + B) +uuj 2 a }, (Ala) 
T Y (u) = D{[u(u + B) + A 2 }(u + C) +uuo\}, (Alb) 
T(u) = D{B-C)uuo Al (Ale) 

where D denotes the function 

B{u) 

^ = {u[u + B{u)} + A 2 }[u + B(u)} + uuj\ ' ( ' 
The extra functions B and C are defined by 

JX.)-<°M3>, CM-ffiffi. (A3) 



(G(»)) 1 
where we have introduced 

G(u) 



(G{u) lR ) ' 



i(u + 7^) + A 2 ] (u + 77?,) + uoj\ ' 



(A4) 



Using that the Laplace transform of f(t)e ±luJAt is given 
by /(u^i'wi), in the case A = it is possible to recuper- 
ate the expressions of Section III- A, Eqs. (|25|l . On the 
other hand, taking u>a = it is straightforward to get 
the results of Section III-C, Eqs. 
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In an effective approximation, Eq. (|18|) . the corre- 
sponding kernels read 



rx(t) 

Ty(t) 

T(t) 

$y(i) 



^){(-) 2 + (— ) 2 cosM}, 
-K'(t) cos[<^t], 

K(t)^{l-co S [<pt}}, 
K(t)—sm[<pt], 



(A5a) 
(A5b) 
(A5c) 

(A5d) 

(A5e) 



where 93 = \Jlo 2 a + A 2 . 

The extra kernels $x (i) and <£>y (i) couples the deriva- 
tive of Sz{t) to the averages Sx(t) and Sy(i) respec- 
tively, i.e., dS z {t)/dt = AS Y (t)- fo dT{$ x (t~T)S x (T) + 
$y(i — t)Sy(t)}. For the exact evolution, these kernels 
vanish. 



APPENDIX B: KERNELS IN THE 
SPIN-BOSON-MODEL NOTATION 



where Q'(t) and Q"{t) are defined by 

Q,{t) = £L 



fc^W coth(^)[l - cos(«,t)], (B7) 



, J(w) . , , 

aw — — su\{wt). 



(B8) 



The exact evolution Eq. I|24ll can be written as in Eq. I|43|) 
with the definitions 



1 A 


(u) = 


)=K ( A a \u) 
K A a \u) = 


= 0, 
0, 


(B9a) 
(B9b) 


K {s) 

^A 


(u) = 


T(u)A 2 [u- 


f r x (u)], 


(B9c) 


V ( a ) 
1 A 


(u) = 


T{u)A[u A 


-T(«)], 


(B9d) 


Y (s) 
1 B 


(u) = 


T{u)A[u + 


Tx(u)}, 


(B9e) 


B 


(u) = 


T{u)A[u A 


-T(u)], 


(B9f) 


Y (s) 
1 B 


(u) = 


T{u)[u + Y 


y(u)], 


(B9g) 


B 


(u) = 


-T(u)[wa 


-T(«)], 


(B9h) 



where we have introduced 



The kernels of the spin-boson model Eqs. (|43|l . in low- 
est order in A read 221 



1 A 


(*) 


~ -AyW(t)Bin[Q"(t)], 


(B6a) 


K {a 

A 


(*) 


~ A 2 Y^(t) S m[Q"(t)l 


(B6b) 


A 


(*) 


~ A 2 ri s) (t)cos[Q"(t)], 


(B6c) 


V (a 
1 A 


(*) 


~ -AY^(t)co S [Q"(t)l 


(B6d) 


K (s 

B 


(*) 


~ Ay«(t), 


(B6e) 


B 


(*) 


~ -Ayj°>(f), 


(B6f) 


v( s 

1 B 


(*) 


~ cos[w j4 ^]e~ Q ' (t) , 


(B6g) 


V ( a 
1 B 


(*) 


~ -sin[o;^t]e- Q 'W, 


(B6h) 



T(u) 



h - T( u )] 2 + [« + r x (w)][u + ry(u)] ' 



(BIO) 



The structure of these kernels is the same as those of 
the spin-boson model in the limit of vanishing A, which 
implies that A only appears through the unitary evolu- 
tion. In fact, in this limit we can approximate Eqs. (|A1() 
by Tx(u) ~ Ty(u) ~ [K(u — iuj A ) + K(u + iuj A )]/2 and 
T(u) ~ [K (u — iu> A ) — K(u + iui A )]/2i. After introduc- 
ing these expressions in Eqs. I|B9|) . it is simple to get 

~ cos[u A t]P (t), and_l^ a) (i) ~ - sm[u A t]P (t). 
Then, disregarding in Eqs. (|B6|I the phase contribution 
proportional to Q"(t), which is valid in a high temper- 
ature limit :29j, a mapping with Eqs. i|B9l) can be done 



after imposing the equality Po(t) = e 



-Q'(t) 
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